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Abstract 

Several methods for solving efficiently the one-dimensional deconvolution prob- 
lem are proposed. The problem is to solve the Volterra equation leu := 
Jq k(t — s)u(s)ds = g(t), < t < T. The data, g(t), are noisy. Of spe- 
cial practical interest is the case when the data are noisy and known at a 
discrete set of times. A general approach to the deconvolution problem is 
proposed: represent k = A{I + S), where a method for a stable inversion of 
A is known, S is a compact operator, and I + S is injective. This method 
is illustrated by examples: smooth kernels k(t), and weakly singular kernels, 
corresponding to Abel-type of integral equations, are considered. A recursive 
estimation scheme for solving deconvolution problem with noisy discrete data 
is justified mathematically, its convergence is proved, and error estimates are 
obtained for the proposed deconvolution method. 



1 Introduction 



In many applications one models the relation between input signal u(t) and output 
signal g{t) by the equation 

ku:= [ k(t- s)u{s)ds:= k*u = g(t), < t < T, (1.1) 



where k(t), given for all t > 0, characterizes the linear system, k*u is the convolution, 
u(t) = k{t) = for t < and the cases T < oo and T = oo are both of interest. 
In practice is measured with some error, so g$(t) is known, \\gs(t) — g(t)\\ < S. 
The norm we use is L 2 (0, oo; e' 2at ) norm, or L 2 (0,T), or L°°(0,T), and the case 
T < oo can be reduced to the case T = oo, as we show below. 

If the operator k in (1.1) is considered as an operator on L°°(0,T), and 
Jq \k(t)\dt < oo, then k is not boundedly invertible, so problem (1.1) is ill-posed. 
One can see this from the formula J* k(t — s)e ms ds — > as n — > oo. 

If T < oo, one sets u(s) = for t > T and defines g(t) for t > T as the left- 
hand side of equation (1.1). If this is done, then (1.1) can be considered as an 
equation on (0, oo) and its solution equals to u(t), the solution of (1.1) on [0, T], 
when t G [0,T]. We assume that k(t) ^ and f Q \k(t)\dt < oo. Then (1.1) has at 
most one solution in L^O, T] if k(t) ^ almost everywhere in [0,T] ([10], p. 327). 
The assumption J °° \k(t)\dt < oo can be replaced without loss of generality by a 
weaker assumption J °°exp(— at)\k(t)\dt < oo, with an arbitrary large fixed a > 0. 
This weaker assumption can be reduced to the original one by changing variables. 

A deconvolution method is a method to construct a stable approximation u$(t) 
of the solution u(t) to equation (1.1), given gg: 

\\u s (t) - u(t) || := T]{8) ^0 as 5^0. (1.2) 

An operator R(S) which constructs such us from g$, u$ = R(S)gs, is called a regu- 
larizer (or a regularizing family, since 5 — > 0) if (1.2) holds. 

There is a large literature on ill-posed problems. General methods for con- 
structing regularizers have been developed. They include variational regularization, 
iterative regularization, method of quasisolutions, etc [9]. In Section 2 the specific 
form of equation (1.1) is used for constructing regularizers for equation (1.1). The 
emphasis is on the causality property of the regularizer. The idea is similar to the 
one in [6] and [3]. In Section 3 a simple general method to construct regularizers 
for equation (1.1) is proposed. This method is practically efficient. It is illustrated 
by two examples in which the results from [2]- [5], and [7] are used. In Section 4 
we investigate a recursive algorithm proposed in [1] for solving equation (1.1) with 
noisy discrete data. Again, the emphasis is on the causality property of the es- 
timate: we use only the data collected up to the time t in order to estimate the 
signal u at this time. If one would use a variational regularization for constructing 
a stable estimation of u(t), one would have to use all the data collected on the full 
time inteval [0, T], and not only on the "current" time interval [0, t\. Our analysis is 
much shorter than in [1] and yields more detailed results. Moreover, In Section 6 we 
discuss briefly a generalization of these results to the case of operator-valued ker- 
nels, which includes, in particular, matrix- valued kernels, that is, systems of Volterra 
equations. In Section 5 proofs are given. 



2 



2 A deconvolution method 



Let K(X) :— J e Xt k(t)dt. By capital letters the Laplace transform is denoted. If 
u solves (1.1) then U(X) = K~ 1 (X)G(X). One has: 

Mt) ■= 7T- / e xt K-\X)G 5 (X)— x -^dX, X = a + ifieC a , (2.1) 

^1 Ja-ioo \Jj + 1J 

where C a is the contour = a, X = a + i/i, m > is a sufficiently large positive 
integer and iV > is a large parameter. We do not show the dependence of us(t) 
on m and N to simplify the notations. We want to prove that one can choose 
N = N(8) -> oo as 8 -> 0, such that (1.2) holds. Note that G 5 (A) = (7(A) + W(A), 
where |W(A)| = |/ °° e~ xt w(t)dt\ < |, the noise iu(i) satisfies the inequality \w(t)\ < 
5, a = 9?A > 0, and G(A) := / °° e~ xt g{t)dt. When one uses L°°(0,T) norm one 
assumes that 



K-\X)G{X) G L\C„), (2.2) 

This is an a priori assumption on u(t). 
We assume throughout this paper that: 

\K(X)\ > c\X\' a , X G C a , a G K. (2.3) 

Here and below c > denote various constants independent of 5 and A/". The 
constant a may be negative, but in many applications a > 0. 



Theorem 2.1 If m > a + 0.5 and (2.3) holds, then there exists N(8) — > oo as 
5^0, swc/i t/iat (1.2) aoWs win L 2 (0, oo; e~ 2<T *) norm. 

If (2.2) and (2.3) hold, and m > a+ 1, then (1.2) holds with L°°(0,T) norm. 

If m > a + 1, (2.3) holds, and 

\K~ 1 (X)G(X)\ < 1 + | C A|1+d , d = const > , A G C a , (2.4) 

then N(S) = O ^ ? q — a + d + 1, rj(S) = 0(6*) if < d < 1, and the norm in 
(1.2) is L°°(0,T) norm. 

Ifd>\, then N(5) = 0(8'^) and n(5) = 0(8^). 

Ifd=l, then N(8) = 0(ex P [M(i + o( i )]); and v{5 ) = 0( _ I ^iL__). 

We sketch proofs in the last Section. 
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3 A general approach to deconvolution 

Suppose the operator k in (1.1) can be decomposed into a sum k := A + B, where 
A~ 1 B := S is compact in the Banach space X, in which k acts, and I+S is boundedly 
invertible, or which is the same by the Fredholm alternative, N(I + S) = {0}, where 
N(A) is the null space of A. In this case / + S is an isomorphism of X onto X, 
11(A) = TZ(k), and 

ku = A(I + S)u = g . (3.1) 
If a regularizer for A is known, then (3.1) can be solved stably by the scheme 

u s = (I + S)- 1 R(5)g s , (3.2) 

and (1.2) holds. 

Since I + S is an isomorphism, the error \\v — vg\\ of the approximation of the 
stable solution of the equation Av = g by the formula v$ = R(5)gs is of the same 
order as \\us — u\\. 

Example 3.1. Let k(t) G C l {Q,T) and assume fc(0) ^ 0. Then, without loss of 
generality, one may assume fc(0) = 1. Write (1.1) as 

ku — u(s)ds + / [k(t — s) — 1] u(s)ds := Au + Bu = g. 
Jo Jo 

Here stable inversion of A is equivalent to stable numerical differentiation of noisy 
data. This problem has been solved in [2] (see also [3] - [5], [7]), and the results of 
these works yield the following theorem: 

Theorem 3.1 Assume that II < m 2 < oo. Then the operator 

g s (t + h(5))-g(t) 



R(S)g 5 



h(5) 

1/2 



( R \ ' 

is a regularizer for the operator A if h(S) = 2 ( . One has: \\R(8)gs — 

u\\rco < 2Jrru5 . 

M ^[0,T-h] — V Z 

In [5] weaker a priori assumption on u is used: ||u|| a < m a , < a < 1, where the 

tt i i ■ i r> i II n \ u (t) ~ U i S )\ I \ I 

Hoelder-space norm is dehned as ||u|| a := sup — ■ ■ V sup \u(t)\ . 

X yiy,x,ye[0,T] \t ~ S\ a 0<t<T 

In this example our method yields the equation 

(I + S)u s = R(S)g s , 

where S is a Volterra operator: Sus := f^k'^t — s)u$(s)ds. Therefore us can be 
easily found by iterations. 
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Example 3.2. Let k(t) = j^^y + m(t), < 7 < l,m(t) G C 1 , 

Au = y * w i Bu = m* u. 

One has A^g = ^ ft jf^. Define 

where R(S) is defined in Theorem 3.1. 

Theorem 3.2 The operator-function Ri(5) is a regularizer for equation Au = g, 
the operator S := A~ X B is compact in L 2 (0,T) and N(I + S) = {0}. 

Compactness of S is clear from its definition. The operator R\(S) is a regularizer for 
A because R(5) is a regularizer for the operator of differentiation. Finally N(I+S) = 
{0} because A~ X B is a Volterra operator. Therefore us can be easily computed by 
our general method and (1.2) holds. 



4 Recursive estimation given discrete noisy data 

Assume that 

gs(nh) := £ n = g(nh) + w n 

are noisy measurements of the data g(t) at the time moments nh, h > is small 
number, \w n \ < S is noise. One wishes to estimate stably u(t), the solution to (1.1), 
given the data £1, ....£„. The following estimation method was essentially proposed 
in [1]. Set v = ^j^ 1 - Define Vj by recursive formulas: 



aVn + / k ( nh ~ s ) dsv i = €n, v = 1 ° . (4.1) 



In this Section we assume that k(t) G L 1 (M + ) fl Lf oc (R + ) and u(t) G Ll c (R + ). Then 
g(t) G Ci oc (R + ), as the following lemma claims. 

Lemma 4.1 Ifu,ke Lf oc (R + ) then k*ue Q oc (R + ). 

Proof. One has 



rt+n 

\(k*u)(t + h) - (k*u)(t)\ < / \k(t + h-s) - k(t - s)\\u\ 

Jo 

/t+h 
\k{t- s)\\u\ds := h + I 2 - 



ds 
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Now, 



/•t+h pt+h 

If< \k(t + h- s)-k(t- s)\ 2 ds / \u\ 2 ds -> as h -> 0, 
Jo Jo 



and 



|jfc(* - s)| 2 ds y |w(s)| 2 ds -> as /i -> 0. 



n Jt 
Lemma 4.1 is proved. □ 

This lemma makes it reasonable to assume that g(t) G Cz oc (IR+). In [1] the case 
g G L 2 oc (M, + ) is discussed, when g is not defined pointwise. It is proposed in [1] to 
use a mollification of g around the points nh instead of using g{nh). However, this 
mollification requires a knowledge of g in a neighborhoods of all points nh, and this 
is an information different from the one assumed at the beginning, namely £i, ....£ n . 
By this reason and because of Lemma 4.1, we assume that g(t) G C/ oc (lR + ), the 
space of functions continuous on any compact subinterval of K+. 
Our assumptions in this section are: 

A) k(t) G L\R + ) n LL(K+) and u(t) G Ll c (R + ), 

B) (2.3) holds, 

C) the union of the spectra of {X(A)}vAec CT does not contain the set {z : z G 
C, 7r — if < arg z < 7r + ip, \z\ < r} , where </? > and r > are arbitrary small fixed 
numbers. 

If the assumptions A), B), and C) hold, then our result is Theorem 4.1 below. 
If, in addition, k and g are Hoelder-continuous, then our result is Theorem 4.2 below, 
which gives the rate of convergence in (1.2). 

We prove that if a = a(5) and h = h(S) are chosen suitably, then the function 
vs(t), defined by the formula: 

v s (t) = Vj for jh<t< (j + l)h, 

approximates stably u(t), so that || ?;,$(£)— u(t)\\ — > as S — > 0, where ||-|| is L°°(0, T)- 
norm. The rate of convergence is estimated in Theorem 4.2 under additional a priori 
assumptions. 

Let t 

au a {t) + I k{t- s)u a {s)ds = g(t), < t < T. (4.2) 
Jo 

The function v${t), defined above, solves the equation 

av s {t)+ [ k{t - s)v s {s)ds = f s {t), (4.3) 
Jo 
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where fs(jh) = for t = jh, and for other values of t the function fs(t) is defined 
as the left-hand side of (4.3): 

f s (t) := g s (nh) - g(t) + / [k(t - s) - k{nh - s)} v s {s)ds + g(t) := tp s (t) + g(t), 

Jo 

where (n — l)h < t < nh, n = 1,2, 3, 

One has 

\Mt)\<S + lg (h) + c(T)\\v 5 \\ lk (h), 

where 

1g(h):=\g(nh)-g(t)\->0, j k {h) := [ \k(t - s) - k(nh - s)\ds -> as h -> 0, 

and (n — < t < nh, n = 1,2, 3 

Denote — u a := w. From (4.2) and (4.3) one gets: w = (a + k)" 1 ^, where 
a + k := al + k, and I is the identity operator. One has ||(q! + k) _1 || < ca~ l , and 
1 < £ + 7 fl (>0 + c(T) 7fc (ft)|M. Therefore: 

K-. a ||<c(T)^±^±^M^M. 

a 

Let us choose /i = h(5) such that 7fc ( fe )+ 7 a( fc ) — > as a — > 0, and a = a(5) so that 
— jt\ ^ as S > . 

a (5) 

Let := u a ($). Then (1.2) holds by Lemma 5.1 below: 

\\us — u\\ = \\u a (s) — u\\ = \\(ce + k) _1 ku — u\\ := e(a) — > as a := a(5) — > 0. 
/,From (4.3) one gets 

= (a + k) _1 + (a + k) _1 ^(n/i) - #(f)] + (a + k) _1 Qu tf , (4.4) 

where t 

Qi^ := / — s) — &(n/i — s)]v$ds. 
Jo 

Since ]| (ct + 1| < % and \\g s (nh) - g(t)\\ < c{T){5 + j g {h)), one has: 

IKa + k)- 1 ^)-^)]! <c{T) 5 + l9{h) ^0 5^0, 

a 

provided that a = a(5) and /i = are chosen so that — > and 7g ^^ — > 
as 5 -> 0. Let a = a(<J). Then ||(a + k^Q^H < ^7 fe (/i)|h||, and if -> as 

5^0, then ||(a + k^QH < -> as 5 -> 0. Thus 



ir»i-i 



i)j= [/-(a + kJ-'Q] 



(o+10 -. 9+o (£±*(!!)) 



7 



where O denotes an element whose norm is O ^ <5+ ^( fe ) j . Therefore, with 

m<5 := (a(5) + k) _1 (7, one gets: 

5 + lg (h) + lk (h) 



\vs -us\\ < c(T)- 



Consequently: 

Ih-ull < \\v S -u s \\ + \\u s -u\\<c(T) S + lfg ^ +lk ^ +e(a)^0 as 5^0. 

a 

(4.5) 

We have proved: 

Theorem 4.1 Assume A), B) 7 and C). Then there exist h = h(8) — > and a = 
a(S) — > snc/i i/iai (4.5) holds. 

If and are Hoelder-continuous, then j g (h) + 7fc(/i) = 0(h b ), < b < 1. 
Put 5 = h b . Then (4.5) can be written as: 

\\v s - u\\ < c\- +e(a)]. (4.6) 
a 

Let us estimate e{ot) = \\u$ — u\\. One has: 

/oo 
|a(a + ^(A))" 1 ! |^- 1 (A)G(A)| d/x 
-oo 



-oo 
r-oo 



/oo 
aKa + i^A))" 1 ! 
-oo 



[l + ^ + ^Vaji+d 

/M pM poo 

+ + :=h + I 2 + h. 
-M J-oo JM 

Using (2.3), one gets: 

h < ca [ M < c«M fl - d . 

x - 7^(1 + 1^)^1 - 

The estimates of I 2 and J 3 are similar. Let us estimate, for example, I 3 : 
Thus 

e(a) < c(aM a ~ d + M~ d ). 

If d > a, then £(«) < ca, Hi^ — u\\ < c[- + £(«)]. Minimizing with respect to a, 
one gets \\v s - u\\ < c5 1/2 if a = S 1/2 . If d < a, then e(a) < c(aM fl_d + M~ d ). 
Minimizing with respect to M, one gets e(a) < ca d / a , Hi^ — it 
Minimizing with respect to a, one gets \ \v$ — u\ \ < c5~ . 
Let us summarize the result: 
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Theorem 4.2 Assume A), B) 7 and C). If g and k are Hoelder- continuous, so that 
lg{h) + Jk(h) = 0(h b ), < b < 1, then \\v$ — u\\ < cS 1 ^ , provided that h = 5^ , 

a i 

a = 0(5 d + a ), and d < a. If d> a, then \ \vs — u\ \ < c5z . 

In [8] a singular perturbation problem was solved for a class of one- and multi- 
dimensional integral equations. The problem we study in Sec. 4 contains a singular 
perturbation problem as a basic component: we are interested in the behavior of 
the operator (a + k) _1 as a — > +0. 



5 Proofs. 

The norm below is L 2 (0, oo; e~ 2at ) norm, it is equivalent to L 2 (0, T) norm on (0, T). 
By the spectrum of a scalar function K(X) we mean the set of its values, and if K(X) 
is an operator-valued function, then its spectrum is defined as usual. 

Lemma 5.1 Let (2.3) hold and assume that the union of the spectra of ^{K(A)}vAec* CT 
does not contain the set {z : z G C, 7r — tp < axgz < 7T + <p, \z\ < r} , where <p > 
and r > are arbitrary small fixed numbers. Then 

e(a) := \\(a + k) _1 kw - u\\ -> as a — > +0 for all u6L 2 (0,oo;e" M ). 



Proof. One has 



{a + k) _1 ku - k^k^ 2 = || (a + k^ak^kuH 2 

= + k) _1 n|| 2 



/oo 
e l ^a(a + K(A))- 1 ^- 1 (A)G'(A)(i/i 
-oo 

i r°° 

= — a 2 |(a + K(A))- 1 | 2 |K- 1 (A)G'(A)| 2 ^:=£(a) ^0 as a^O, 
A = a + i/i G C CT . 

Here we have used: 1) Parseval's equality; 2) the asumption u G L 2 (0, oo; e~ 2<T * 
which is equivalent to K~ 1 (X)G(X) G L 2 (C a ); 3) the estimate sup a(a + K(X))^ 1 

< c which follows from the assumption about the range of K(\) on C a \ 4) estimate 
(2.3); and 5) the dominated convergence theorem. Lemma 5.1 is proved □ 

Lemma 5.2 Under the assumptions of Lemma 5.1 one has: 



c5 

(a + k)^ 1 ^ — u\\ < h s(a), e(ot) ^ as a — ■> 0, c = const > . 

a 
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Proof. One has 

IKa + k)" 1 ^-^ < || (« + k) - x (^ - ^) || + || (« + k)" 1 ^ - u|| < IKa + k)- 1 )!!^^). 

By Lemma 5.1, e(a) — > as a — > 0, and ||(a + k) _1 || < caT 1 , where c =const> 
depends on tp, as follows from the proof of Lemma 5.1 and from the estimate 



sup |(a + K(A))" 1 | < 



\ec a 



a 



Lemma 5.2 is proved. □ 

Proof of Theorem 2.1. If u$ is defined in (2.1) and 
then Parseval's equality yields 



is L 2 (0,oo;e 2at ) norm, 



\us — u\ 
1 



1 

27 



K-\X)Gs(\) 



< 



7T 



^-'(AJGfA)! 1 



(* + 1)' 



djji + 



^ _1 (A)G(A) 
o 2 



d\x 



2m 



JV 



+ 1 



/l + $ 2 h 



where we have used the formulas Gs = G + W, \W\ < -. 

If N — > oo then ix = h{N) — > by the dominated convergence theorem. Let us 
estimate I2 assuming (2.3) and taking m — a > 0.5: 



I 2 < cN 

< c 



2m 



{a 2 + fi 2 ) a dfi 



< cN 2m 



dfx 



(N + 



-00 [(a + N) 2 + /i 2 ]- " " 7^ [(a + N) 2 + /i 
^ r ± < cA r- + i 

(T )2m-2a-l J ^ (J + ^m-a - 



2]m—a 



Thus, if 2a + 1 > 0, then, using the estimate (x + y) 1 / 2 < x 1 ! 2 + y 1 / 2 , rr, y > 0, one 
gets: 

\\u s - u\\ < 5c 1/2 N a+1 z + I*(N) := 7](S, N) . 

Minimizing r](5, N) with respect to iV for a fixed 8, denoting the minimizer by N(S), 
N(S) — > 00, as 5 — > 0, and the minimum by 77(0") := r}(8, N(5)), one gets 77(0") — > 
as 5 — > 0. Thus, (1.2) is proved with L 2 (0,oo; e~ 2 °"*) norm. □ 

If L°°(0,T) norm is used for \\u(t) — v>s(t)\\, then one gets (1.2) if assumptions 
(2.2) and (2.3) are used. Namely, 



sup \u$(t) — u(t)\ < 

0<t<T 



exp(o\T) 
2^ 



\K- l G s (l + ±) 



< c(T)(Ji + J 2 ) := 7/, 
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where 



and 



/oo \ 
\K- 1 G\\(l + -)- m -\\dp, 
-oo ^ 

J 2 ;= J 2 (N) = - / + 

If K~ X G G L 1 (C cr ), then liniTv^oo <A = by the dominated convergence theorem. If 
m > a + 1 and (2.3) holds, then 

J 2 < cSN a+1 . 

Thus (1.2) holds with L°°(0, T)-norm provided that K^G E L 1 ^), m > a + 1 
and (2.3) holds. 

If (2.4) holds, one can get a rate of decay. Namely 

pM p 

Ji= I + I =3i+32, 

where 



M J\p\>M 



/M \ 
^(i+iAr i )- i Ki+-)-™-i|dii 

■ J^ (1 + lAI^ 1 )" 1 ^^ < cM^N' 1 , if < d < 1, 



< C 



Ji < ciV" 1 if d > 1, 
Ji < c ^f~ if rf= 1, 



and 



32 <c (l + |Ar i )- 1 |(l + ^)— -l\dpt<c / (1 + IAO- 1 < cM- d , 

c > stands for various constants, and A = a + Thus 

Ji < c(M l - d N~ l + M~ d ) if < d < 1, 

Ji < c(iV- 1 + M~ d ) if d > 1, 

J^c^ + M^) ifd=l. 

If < d < 1 then choose M = N and get Ji < cA^ d . Therefore, if m > a + 1 and 
(2.3) and (2.4) hold, then 

Ji + J2 < c(5N a+1 + N- d ) if < d < 1. 
11 



Minimizing with respect to N for a fixed S > 0, one gets the minimizer iV = N(5) = 
0(5~ T +^+ 1 ) and the estimate i] < 0(5 T +^+ 1 ) if < d < 1. If d > 1 then choose 
M d = N and get Ji < cN" 1 , Ji + J 2 < c(SN a+1 + A^ 1 ), iV(<J) = 0(<T^), and 
the estimate 77 < 0(J2+^). If d = 1 then choose M = N and get J\ < c 1 - 2 ^- 



Jx + J 2 < c(SN a+1 + N(S) 



0(exp[M ( 



N 



Ijogjl 



a+2 



1 + o(l)]), and the estimate rj < 



cxp 



I log 5 
2+a 



(1+°(1))] ' 



Theorem 2.1 is proved. □ 



Remark 5.3. If k(t) > and k(t) 
If 3?K(A) > then a(a + K(X)Y 
assumption of Lemma 5.1 holds with ip = | and r = 00. 



monotonically, then KfT(A) > for 3?A > 0. 
< 1. Condition 3WT(A) > implies that the 



6 Generalizations. 

Most of our results and proofs remain valid for operator- valued functions k(t), in 
particular for matrix- valued kernels, that is, for systems of Volterra equations. Let 
k{t) be an operator in a Banach space, and K(X) be its Laplace transform. If 
one replaces the absolute values by the norms in (2.3), (2.4) and elsewhere in the 
proofs, then one gets Theorems 2.1, 3.1, 3.2, and 4.1 and lemmas 5.1 and 5.2 with 
operator- valued k(t). 

Acknowledgement. AGR thanks Prof. L. Pandolfi and Dr. F. Fagnani for 
discussions. Their paper [1] was useful for the results in Sec. 4 of this paper. 
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